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In order to find a way to have a better formulation for numerical evolution of the Einstein equa- 
tions, we study the propagation equations of the constraints based on the Arnowitt-Deser-Misner 
formulation. By adjusting constraint terms in the evolution equations, we try to construct an 
"asymptotically constrained system" which is expected to be robust against violation of the con- 
straints, and to enable a long-term stable and accurate numerical simulation. We first provide 
useful expressions for analyzing constraint propagation in a general spacetime, then apply it to 
Schwarzschild spacetime. We search when and where the negative real or non-zero imaginary eigen- 
values of the homogenized constraint propagation matrix appear, and how they depend on the choice 
of coordinate system and adjustments. Our analysis includes the proposal of Detweiler (1987), which 
is still the best one according to our conjecture but has a growing mode of error near the horizon. 
Some examples are snapshots of a maximally sliced Schwarzschild black hole. The predictions here 
may help the community to make further improvements. 



I. INTRODUCTION 



Computing the Einstein equations numerically is a necessary direction for general relativity research. This approach, 
so-called numerical relativity, already gives us much feedback to help understand the nature of strong gravity, such 
as singularity formations, critical behavior of gravitational collapses, cosmology, and so on. Numerical relativity will 
definitely be required in unveiling gravitational wave phenomena, highly relativistic astrophysical objects, and might 
be the only tool to discuss different gravitational theories or models [1]. 

There are several different approaches to simulate the Einstein equations, among them the most robust is to apply 
3-1-1 (space -|- time) decomposition of spacetime, which was first formulated by Arnowitt, Deser and Misner (ADM) 
[2] (we refer to this as the "original ADM" system). A current important issue in the 3-1-1 approach is to control 
long-term stability of time integrations in simulating black hole or neutron star binary coalescences. So far, the use of 
fundamental variables (5^ , -ftTy ), the internal 3-metric and extrinsic curvature, has been the most popular approach 
[3] (we refer to this as the "standard ADM" system). In recent years, however, many groups report that another re- 
formulation of the Einstein equations provides more stable and accurate simulations. We try here to understand these 
efforts through "adjusting" procedures to the evolution equations (that we explain later) using eigenvalue analysis of 
the constraint propagation equations. 

At this moment, there may be three major directions to obtain longer time evolutions. The first is to use a 
modification of the ADM system that was developed by Shibata and Nakamura [4] (and later remodificd by Baumgarte 
and Shapiro [5]). This is a combination of the introduction of new variables, conformal decompositions, rescaling the 
conformal factor, and the replacement of terms in the evolution equation using momentum constraints. Although 
there are many studies to show why this re-formulation is better than the standard ADM, as far as we know, there is 
no definite conclusion yet. One observation [6] pointed out that to replace terms using momentum constraints appears 
to change the stability properties. 

The second direction is to re-formulate the Einstein equations in a first-order hyperbolic form [7] [33]. This 
is motivated from the expectation that the symmetric hyperbolic system has well-posed properties in its Cauchy 
treatment in many systems and also that the boundary treatment can be improved if we know the characteristic speed 
of the system. We note that, in constructing hyperbolic systems, the essential procedures are to adjust equations 
using constraints and to introduce new variables, normally the spatially derivatived metric. Several groups report 
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that the hypcrbohc formulation actually has advantages over the direct use of ADM formulation [9, 10, 11, 12]. 
However it is also reported that there arc no drastic changes in the evolution properties between hyperbolic systems 
(weakly/strongly and symmetric hyperbolicity) by systematic numerical studies by Hern [13] based on Prittelli-Reula 
formulation [14], and by the authors [15] based on Ashtekar's formulation [16, 17, 18] [34]. Therefore we may say 
that the mathematical notion of "hyperbolicity" is not always applicable for predicting the stability of numerical 
integration of the Einstein equations (See also §4.4.2 in [1]). It will be useful if we have an alternative procedure to 
predict stability including the effects of the non-principal parts of the equation, which are neglected in the discussion 
of hyperbolicity. 

The third is to construct a robust system against the violation of the constraints, such that the constraint surface 
is the attractor. The idea was first proposed as "A-system" by Brodbeck et al [21] in which they introduce artificial 
flow to the constraint surface using a new variable based on the symmetric hyperbolic system. This idea was tested 
and confirmed to work as expected in some test simulations by the authors [22] (based on the formulation developed 
by [23]). Although it is questionable whether the recovered solution is true evolution or not [24], we think that 
to enforce the decay of errors in its initial perturbative stage is the key to the next improvements. Actually, by 
studying the evolution equations of the constraints (hereafter we call them constraint propagation) and by evaluating 
eigenvalues (amplification factors, AFs) of constraint propagation in its homogenized form, we found that a similar 
"asymptotically constrained system" can be obtained by simply adjusting constraints to the evolution equations, even 
for the ADM equations [25]. 

The purpose of this article is to extend our previous study [25] in more general expressions and also to apply 
the systems to spacetimc which has non trivial curvature, Schwarzschild black hole spacetime. The actual numerical 
simulations require many ingredients to be considered such as the choice of integration schemes, boundary treatments, 
grid structures and so on. However, we think that our approach to the stability problem through an implementation 
to the equations is definitely one of the aspects that should be improved. 

Adjusting evolution equations is not a new idea. Actually the standard ADM system for numerical relativists [3] is 
adjusted from the original one [2] using the Hamiltonian constraint (See Frittelli's analysis on constraint propagation 
between the original and standard ADM formulations [26]). Detweiler [27] proposed a system using adjustments so 
that the L2 norm of constraints may not blow up (we will study this later in detail). Several numerical relativity 
groups recently report the advantages of the adjusting procedure with a successful example [10, 11]. We try here to 
understand the background mathematical features systematically by using AFs of constraint propagation. 

We conjecture that we can construct an asymptotically constrained system by evaluating the AFs in advance. 
Especially, we propose that if the amplification factors are negative or pure-imaginary, then the system becomes more 
stable. This conjecture is already approved by numerical studies [22] in the Maxwell equations and in the Ashtekar 
version of the Einstein equations on the Minkowskii background. In this article, we show when and where such a 
feature is available using the ADM equations for Schwarzschild spacetime. Although the analysis here is restricted 
to the fixed background, we believe that this study clarifies the properties of adjustments and indicates further 
possibilities. 

In section II, we describe our idea of "adjusted systems" generally, and their application to the ADM formulation. 
We then explain how we calculate the "amplification factors" in the spherically symmetric spacetime in §111, and show 
several examples in §IV. The concluding remarks are in §V. Appendix A summarizes the constraint propagation 
equations of the ADM formulation generally, which may be useful for further studies. Appendix B explains briefly a 
way to obtain maximally sliced Schwarzschild spacetime that is used to obtain Fig. 7. 



II. ADJUSTED SYSTEMS 



A. Procedure and background — general discussion — 

We begin with an overview of the adjusting procedure and the idea of background structure, which were described 
in our previous work [22, 25]. 

Suppose we have a set of dynamical variables u°'{x'^,t), and their evolution equations 

atu« = /(««,ai««,---), (2.1) 

and the (first class) constraints 

C«(w°,ai«»,---)«0. (2.2) 

For monitoring the violation of constraints, we propose to investigate the evolution equation of (constraint prop- 
agation) , 

5tC" =5(C",5iC",---). (2.3) 
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(We do not mean to integrate (2.3) numerically, but to evaluate them analytically in advance.) Then, there may be 
two major analyses of (2.3); (a) the hyperbolicity of (2.3) if (2.3) forms a first order form, and (b) the eigenvalue 
analysis of the whole RHS in (2.3) after certain adequate homogenized procedures. 

If the evolution equations form a first-order system, the standard hyperbolic PDE analysis is applicable. The 
analysis is mainly to identify the level of hyperbolicity and to calculate the characteristic speed of the system, from 
eigenvalues of the principal matrix. We expect mathematically rigorous well-posed features for strongly hyperbolic 
or symmetric hyperbolic systems. Also we know the characteristic speeds suggest the satisfactory criteria for stable 
evolutions; such as when they are real, under the propagation speed of the original variables, u", and/or within the 
causal region of the numerical integration scheme that is applied. 

For example, the evolution equations of the standard/original ADM formulation, (2.1) [(2.7) and (2.8)], do not form 
a first-order system, while their constraint propagation equations, (2.3) [(2.11) and (2.12)], form a first-order system. 
Therefore, we can apply the classification on the hyperbolicity (weakly, strongly or symmetric) to the constraint 
propagation equations. However, if one adjusts ADM equations with constraints, then this first-order characters will 
not be guaranteed. Another problem in the hyperbolic analysis is that it only discusses the principal part of the 
system, and ignores the rest. If there is a method to characterize the non-principal part, then that will help to clarify 
our understanding of stability. (This importance was mentioned in [1].) 

The second analysis, the eigenvalue analysis of the whole RHS in (2.3), may compensate for the above problems. 
We propose to homogenize (2.3) by a Fourier transformation, e.g. 

dtC" = giC") = M^ffCf^, 

where C{x,ty = j C{k,ty exp{ik ■ x)(fk, (2.4) 

then to analyze the set of eigenvalues, say As, of the coefficient matrix, M"^, in (2.4). We call As the amplification 
factors (AFs) of (2.3). As we have proposed and confirmed in [22]: 

Conjecture: 

(a) If the amplification factors have a negative real-part (the constraints are forced to be diminished), 
then we see more stable evolutions than a system which has positive amplification factors. 

(b) If the amplification factors have a non-zero imaginary-part (the constraints are propagating away), 
then we see more stable evolutions than a system which has zero amplification factors. 

We found heuristically that the system becomes more stable when more As satisfy the above criteria [15, 22] [35]. We 
remark that this eigenvalue analysis requires the fixing of a particular background spacetime, since the AFs depend 
on the dynamical variables, w°. 

The above features of the constraint propagation, (2.3), will differ when we modify the original evolution equations. 
Suppose we add (adjust) the evolution equations using constraints 

= /(«", diu'', ■■■)+ F{C", d,C", ■ ■ •), (2.5) 

then (2.3) will also be modified as 

dtC" = g{C", diC,- ■ ■) + G{C",diC'', ■ ■ ■). (2.6) 

Therefore, the problem is how to adjust the evolution equations so that their constraint propagations satisfy the above 
criteria as much as possible. 

B. Standard ADM system and its constraint propagation 

We start by analyzing the standard ADM system. By "standard ADM" we mean here the most widely adopted 
system, due to York [3], with the evolution equations, 

^a^J = -2aif„-+V,/3, +V,A, (2.7) 
dtK,, = aB^ +aKK,,-2aK,kK^,-Viyja 

+ {ViP'')Kkj + {VjP'')Kki + P'VkKij, (2.8) 
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and the constraint equations, 

H := R^^^ + - KijK'\ (2.9) 
Mi := VjK^i-ViK, (2.10) 

where (7ij, Kij) arc the induced three-metric and the extrinsic curvature, (a, (3i) are the lapse function and the shift 

(3) 

covector, Vj is the covariant derivative adapted to 7^, and is the three- Ricci tensor. 

The constraint propagation equations, which are the time evolution equations of the Hamiltonian constraint (2.9) 
and the momentum constraints (2.10), can be written as 

^^n = f3^{djH) + 2aKn-2aY^{diMj) 

+a(5amfc)(27"'7'^' - 7'"'7'^)X, - ^Y^{dja)Mi, (2.11) 
dtMi = -{l/2)a{d,n)-{d,a)n + (3'{djM,) 

+aKM^ - l3''-f'\daik)M, + [dAh^'M,. (2.12) 

Further expressions of these constraint propagations are introduced in Appendix A of this article. 

C. Adjustment to ADM evolution equations and its effects on constraint propagations 

Generally, we can write the adjustment terms to (2.7) and (2.8) using (2.9) and (2.10) by the following combinations 
(using up to the first derivatives of constraints for simplicity in order to include Detweiler's case, see the next 
subsection), 

adjustment term of daij : +PijH + Q^ijMk + P^^j {"^kH) + q^\j {VkMi), (2.13) 
adjustment term of dtKij : +RijH + S^ijMk + r^ij {VkH) + s^\j {VkMi), (2.14) 

where P, Q, R, S and p, q, r, s are multipliers (please do not confuse Rij with three Ricci curvature that we write as 

According to this adjustment, the constraint propagation equations are also modified as 

dtH = (2.n) + i7{""(2.13) + i?r"5*(2-13) + i?3^""a,aj (2.13) + Hf "(2.14), (2.15) 
dtM^ = (2.12) + Mii""(2.13) + Af2j^''""aj(2.13) + M3/""(2.14)+M4,J""aj(2.14). (2.16) 

with appropriate changes in indices. (See detail in Appendix A. The definitions of • • • , Mi, ■ ■ ■ are also there.) 

D. Examples of adjustments 

We show several examples of adjustments here. 

The first test is to show the differences between the standard ADM [3] and the original ADM system [2]. In terms 
of (2.13) and (2.14), the adjustment, 

Rij = Kpa-yij, (2.17) 

will distinguish two, where Kp is a constant and set the other multipliers zero. Here Kp = corresponds to the standard 
ADM (no adjustment), and Kp = —1/4 to the original ADM (without any adjustment to the canonical formulation 
by ADM). As one can check the adjusted constraint propagation equations in Appendix A, adding Rij term keeps 
the constraint propagation in a first-order form. Prittelli [26] (see also [25]) pointed out that the hyperbolicity of 
constraint propagation equations is better in the standard ADM system. 

The second example is the one proposed by Detweiler [27]. He found that with a particular combination (see below), 
the evolution of the energy norm of the constraints, Ti.^ + A4^, can be negative definite when we apply the maximal 
slicing condition, K = 0. His adjustment can be written in our notation in (2.13) and (2.14), as 



P,j = -KLa^7«j, (2.18) 

Rij = nLa\K,j - (l/3)X7,j), (2.19) 

S\j = tiLa^[i{d(,a)5';) - {dia)^,,-i''% (2.20) 

s'-'ij = «La^[<5(V;.)-(l/3)7.,7''], (2-21) 
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everything else is zero, where kl is a constant. Dctweilcr's adjustment, (2.18)-(2.21), does not put constraint propaga- 
tion equation to a first order form, so we cannot discuss hyperbolicity or the characteristic speed of the constraints. We 
confirmed numerically, using perturbation on Minkowskii spacetimc, that Detweiler's system provides better accuracy 
than the standard ADM, but only for small positive kl- See the Appendix of [25]. 

There are infinite ways of adjusting equations. Actually one of our criteria, the negative real AFs, requires breaking 
the time-symmetric features of the standard ADM system (when we apply them to the time-symmetric background 
metric). One observation by us [25] is that such AFs are available if we adjust the terms which break the time reversal 
symmetry of the evolution equations. That is, by reversing the time {dt —dt), there are variables which change 
their signatures [Kij,dt"fij,A4i,- ■ ■], while not [gij,dtKij,7i,- ■ ■]. (Here wc arc assuming the 3-metric 7^ has plus 
parity for time reversal symmetry). So that the multiplier, for example, Pij in (2.13) is desired to have plus parity in 
order to break the minus parity of dt'fij equation. 

We list several combinations of adjustments in Table I and Table III. (Table I shows the ones wc have plotted, 
and the list of plots is in Table II). We defined the adjustment terms so that their positive multiplier parameter, 
«; > 0, makes the system better in stability according to our conjecture. (Here better means in accordance with our 
conjecture in §11 A). In this article, we do not discuss the ranges of the effective multiplier parameter, k, since the 
range depends on the characteristic speeds of the models and numerical integration schemes as we observed in [22] . 

We show how these adjustments change the AFs in Schwarzschild spacetime in §IV. 

III. CONSTRAINT PROPAGATIONS IN SPHERICALLY SYMMETRIC SPACETIME 

From here, we restrict our discussion to spherically symmetric spacetime. We introduce the violation of constraints 
as its perturbation using harmonics. According to our motivation, the actual procedure to analyze the adjustments is to 
substitute the perturbed metric to the (adjusted) evolution equations first and to evaluate the according perturbative 
errors in the (adjusted) constraint propagation equations. However, for the simplicity, we apply the perturbation to 
the pair of constraints directly and analyze the effects of adjustments in its propagation equations. The latter, we 
think, presents the feature of constraint propagation more clearly for our purposes. 

A. The procedure 

The discussion becomes clear if we expand the constraint := {H,A4i)^ using vector harmonics, 

= ^ {A'^'aim + B'^hm + C'^Clm + D^^^dlJ) , (3.1) 



where we choose the basis as 



a/™(0,(^) ^ (F/„, 0,0,0)^, (3.2) 
bi^{e,^) = {Q,Yi^,Q,Qf, (3.3) 

cim{e,'fi) = -j=^==%0,deYir,„d^Yi^f , (3.4) 
1/ lyl + 1) 

di^{e,^) = —I=={Q,Q,-^d^Yi^,smedeYimf, (3.5) 
y«(< + l) smfy 

and the coefficients A'™, • • • , are functions of {t, r). Here Yi^ is the spherical harmonic function, 



YUe,^) = (-l)('»+NI)/y + ^) || ^ p-(cosg)e-^. (3.6) 
The basis (3.2)-(3.5) are normalized so that they satisfy 

{c^,c,) = j^ d<pj^ ciCpT^'^psmede, (3.7) 

where rj'^'' is Minkowskii metric and the asterisk denotes the complex conjugate. Therefore 

A'" = (a'(^),a), dtA'"' = {o}^ydtC,), etc. (3.8) 



6 



In order to analyze the radial dependences, we also express these evolution equations using the Fourier expansion 
on the radial coordinate, 

^Im ^ ^Jm ^^^^ gifer ^g g^ 

k 

So that we can obtain the RHS of the evolution equations for • • • , in a homogeneous form. 

B. Expression for the standard ADM formulation 

We write the spherically symmetric spacetime using a metric, 

+-frrit,r)dr^ +'yee{t,r){de^ +sm'^0dip'^), (3.10) 

where a and Pr are the lapse function and the shift vector, jrr and ^0$ are also interpreted as 3-metric in 3+1 
decomposition, so that the relevant extrinsic curvature, Kij, become Kij = dieLg{Krr, K00, K00 sin^ 9) and its trace 

becomes K = Krrhrr + '^Keg/-^00. 

According to the procedure in the previous section, we obtain the constraint propagation equations for the standard 
ADM formulation, (2.11) and (2.12), in the following form: 

dtA^"'{t,r) = 2aKA^"'+p''{drA^"') 

-2oL^^\^l{l + 1)C'™, (3.11) 
dtB^"'{t,r) = -{l/2)adrA^'" ^ {dra)A'"' 

+ {aK - /3^7''(9r7rr) + {drPr)YlB^"' + /^''^^S''", (3.12) 

dtC^"'(t,r) = _ ^^^(^ + ^) j^lm ^ g^J^C^m ^f^rA(jlm_^ drC^""), (3.13) 

2r r 

dtD^'^{t,r) = aifZ?'™ + /3''(il?'™ + a,!?'"), (3.14) 



and then 



|i'(^)(i) = {2aK + iW)A'i;:) 

+ ( - 2ikaY'' + aidr^rrWY'^ - 2a{dr^00)^'' - AY''{dra))B\t) 
-2a^''r^M+^)C\Tv (3-15) 

JtB[k){t) = {-i^k/2)a - idra))A% + {aK - Y'^ (drjrr) + {dr^W^^ + iknB[^), (3.16) 

lc\t){t) = -^^^^f^^w + {aK+^+ iW)C\T), (3.17) 
jD[u){t) = {aK+^+iW)D\i:y (3.18) 

There is no dependence on m. We see that the expressions are equivalent to the case of flat background spacetime 
[25] when we take I = and r — > oo. Therefore our results also show the behavior of the flat background limit in its 
large r limit. 

C. Example: original ADM formulation 

We only present here one example, the comparison between the standard and original ADM systems. By substituting 
(2.17) into (2.15) and (2.16), we obtain 

dtU = {2AI) + AniaKT-L, (3.19) 
dtMi = (2.12) - 2Kia{din) - 2Ki{dia)n. (3.20) 
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In the spherically symmetric spacetime, this adjustment affects (3.11)-(3.14) as follows: 

9tA''"(t,r) = (3.11) +4Kia'XA'™, 



dtC^'^{t,r) = (3.13) 



After homogenization, (3.15)-(3.18) become 



(3.14). 



(3.21) 
(3.22) 

(3.23) 
(3.24) 



di 



\ ^(fe) / 




f AmaK 0^ 

-2K.iaik - 2K,i{dra) 

-2Kia^l{l + l)/r 

V 0, 



(3.25) 



The eigenvalues (AFs) A' of the RHS matrix of (3.25) can be calculated by fixing the metric components including 
the gauge. For example, on the standard Schwarzschild metric (4.1), they are 



A' = (0,0,7^,-7^) 



a = -k^ + 



4Mk^r^{r - M) + 2M(2r - M) + l{l + l)r(r - 2M) + ikr{2r^ - 3Mr - 2M^) 



for the choice of ki = (the standard ADM), while they are 

M(2r - M) + irkM{2M - r) 



A' = (0,0, A-Vfe), 
for the choice of ki = — 1/4 (the original ADM). 



b = 



(3.26) 



(3.27) 



D. Our analytic approach 

The above example is the simplest one. In the next section, we will show the AFs of the adjusted systems shown 
in Table I. We found that to write down the analytical expressions of them is not a good idea due to their length. 
We will therefore plot AFs to see if the real parts become negative, or if the imaginary parts become non-zero or not. 

We used the symbolic calculation software, Mathematica and Maple, and made plots by checking two independent 
outputs. These scripts are available upon request. 



IV. CONSTRAINT PROPAGATIONS IN SCHWARZSCHILD SPACETIME 



A. Coordinates 



We present our analysis of the constraint propagation equations in Schwarzschild black hole spacetime, 

2M dr^ 

ds^ = —(1 )dt^ + - — — h r^dfl^, (the standard expression) (4.1) 

where M is the mass of a black hole. For numerical relativists, evolving a single black hole is the essential test 
problem, though it is a trivial at first sight. The standard expression, (4.1). has a coordinate singularity at r = 2M, 
so that we need to move another coordinate for actual numerical time integrations. 
One alternative is the isotropic coordinate, 

= -^TXITT^)'^*' + (1 + ^)^ldrl, + rljn% (the isotropic expression) (4.2) 

which is given by the coordinate transformation, r = (1 + M /2riso)'^Tiso- Here riso = M/2 indicates the minimum 

throat radius of the Einstein-Rosen bridge. Bernstein, Hobill and Smarr [28] showed a systematic comparison for 
numerical integration schemes by applying the coordinate transformation, riso = Me^""^ /2, further to (4.2). 
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The expression of the ingoing Eddington-Finkelstein (iEF) coordinate has become popular in numerical relativity 
for treating black hole boundaries as an excision, since iEF penetrates the horizon without a coordinate singularity. 
The expression is, 

ds^ = -(1 —)dtfEF + —dtiEpdr + (1 + —)dr^ + r^d^^ (the iEF expression) (4.3) 

which is given by tiEP = t + 2Aflog(r — 2M) and the radial coordinate is common to (4.1). The geometrical 
interpretation of the iEF coordinate system is that in addition to having a timelike killing vector, the combination of 
timelike and radial tangent vectors dt — Or remains null. 

Another expression we test is the Painleve-GuUstrand (PG) coordinates. 



ds^ = -\^ —jdrpG + 2^—dtpGdr + dr^+r^dQ.% (the PG expression) (4.4) 

which is given by tpc = t+ VSMr — 2M log{(-\/r/2M + l)/(-\/r/2M — 1)} and the radial coordinate is common to 
(4.1). The PG coordinate system can be viewed as that anchored to a family of freely moving observers (time-like) 
starting at infinity with vanishing velocity [29]. 

The latter two (iEF/PG) are also different from the former (standard/isotropic) two in the point that their extrinsic 
curvature on the initial slice {tiEF = tpG = 0) is not zero. To conclude first, the effects of adjustments are similar 
both between the former and between the latter. 

In Table. II, we list the combinations (adjustments and coordinates) we plotted. 



B. in the standard SchwEirzschild coordinate 



We show first the case of the standard Schwarzschild coordinate, (4.1), since this example provides a basic overview 
of our analysis. The cases of the isotropic coordinate, (4.2), will be shown to be quite similar. 

In Fig. 1(a), the amplification factors (AFs, the eigenvalues of homogenized constraint propagation equations) of 
the standard ADM formulation are plotted. The solid lines and dotted lines with circles are real parts and imaginary 
parts of AFs, respectively. (The figure style is common throughout the article.) They arc four lines each, but as we 
showed in (3.26), two of them are zero. The plotting range is 2 < r < 20 in Schwarzschild radial coordinate. The AFs 
at r = 2 are ±^/3/8 and 0. The existence of this positive real AF near the horizon is a new result which was not seen 
in the flat background [25]. Wc show only the cases with / = 2 and k = 1, because we judged that the plots of Z = 
and other A;s are qualitatively the same. 

The adjustment (2.17) with kf = —1/4 returns the system back to the original ADM. AFs are (3.27) and we plot 
them in Fig. 1(b). We can see that the imaginary parts are apparently different from those of the standard ADM 
[Fig. 1(a)]. This is the same feature as in the case of the flat background [25]. According to our conjecture, the 
non-zero imaginary values are better than zeros, so we expect that the standard ADM has a better stability than the 
original ADM system. Negative kf makes the asymptotical real values finite. If we change the signature oi kf, then 
AFs are as in Fig. 1(c). The imaginary parts become larger than kf = 0, that indicates the constraint propagation 
involves higher frequency modes. 

The adjustment proposed by Dctweiler, (2.18)-(2.21), makes the feature completely different. Fig. 2(a) and (b) are 
the case of Kl = 1/4 and 1/2 and (c) is of the different signature of k^. A great improvement can be seen in both 
positive kl cases where all real parts become negative in large r. Moreover all imaginary parts are apart from zero. 
These are the desired features according to our conjecture. Therefore we expect the Detweiler adjustment has good 
stability properties except near the black hole region. The AF near the horizon has a positive real component. This is 
not contradictory with the Detweiler's original idea. The idea came from suppressing the total L2 norm of constraints 
on the spatial slice, while our plot indicates the existence of a local violation mode. The change of signature of kl 
can be understood just by changing the signature of AFs, and this fact can also be seen to the other examples. 

We next show that the partial adjustments of Detweiler's are also not so bad. Fig. 3 (a), (b) and (c) are the cases 
of adjustment that are of only (2.18), (2.20) and (2.21), respectively. [The contribution of (2.19) is absent since 
Kij = in the Schwarzschild coordinate.] By comparing them with Fig. 1(a) and Fig. 2(b), wc see the negative real 
parts in large r are originated by (2.18) and (2.21), while the former two real parts remain zero. Fig. 3(d) used only 
Pij = —KLCe'fij and everything else is zero, which is a minor modification from (2.18). The contribution is similar, 
and can be said to be effective. 
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C. in isotropic/iEF/PG coordinates 

We next compare AFs between different coordinate expressions. The first test is for the standard ADM formulation. 
Fig. 4 shows AFs for (a) the isotropic coordinate (4.2), (b) the iEF coordinate (4.3), and (c) tlie PG coordinate (4.4). 
All plots are on the time slice of < = in each coordinate expression. [See Fig. 1(a) for the standard Schwarzschild 
coordinate.] We see that Fig. 1(a) and Fig.4(a) are quite similar, while Fig. 4(b) and (c) are qualitatively different 
from the former. This is because the latter expressions (iEF/PG) are asymmetric according to time, i.e. they have 
non-zero extrinsic curvature. 

We note that the constraint propagation equations are invariant for spatial coordinate transformation, but AFs are 
invariant only for linear transformation. This explains the differences between Fig. 1(a) and Fig. 4(a), although these 
are not significant. 

We notice that while some AFs in iEF/PG remain positive [Fig. 4(b) and (c)] in large r region, that their nature 
changes due to the adjustments. Fig. 5 is for the Detweiler-type adjustment, (2.18)-(2.21), for the isotropic/iEF/PG 
coordinate cases. [See Fig. 2(b) for the standard Schwarzschild coordinate.] Interestingly, all plots indicate that all real 
parts of AFs arc negative, and imaginary parts are non-zero (again except near the black hole region). By arranging 
the multiplier parameter, for the iEF/PG coordinates, there is a chance to get all negative real AFs outside the black 
hole horizon. For example, for iEF (PG) coordinate all the real-part goes negative outside the black hole horizon if 
Kl > 3.1(1.6), while large kl may introduce another instability problem [22]. 

Fig. 6 shows the adjustment No. 4 in Table. I, which was used in the test of PennState group [11]. The main 
difference from above is that the adjustment here is only for the radial component of the extrinsic curvature, Krr- 
The numerical experiments in [11] show better stability for positive Kfj,, a fact that can be seen also from Fig. 6: 
positive Kfj, produces negative real AFs. [See Fig. 4(b) (c) for the standard ADM case.] 

Such kinds of test can be done with other combinations. In Table. Ill, we listed our results for more examples. 
The table includes the above results and is intended to extract the contributions of each term in (2.13) and (2.14). 
The effects of adjustments (of each k > case) to AFs are commented upon for each coordinate system and for 
real/imaginary parts of AFs, respectively. These judgements are made at the r ~ O(IOM) region on their t = slice. 
We hope this table will help further numerical improvements for the community. 

D. in mcLximally-sliced evolving Schwarzschild spacetime 

So far, our discussion is limited to one 3-hypersurface. Generally speaking, such an initial-value like analysis may 
not be enough to determine what combination of adjustments, coordinate system, and gauge conditions are suitable 
for the numerical evolution problem. Here as the first further step, we show our analysis of several snapshots of 
maximally sliced Schwarzschild spacetime. 

The so-called maximal slicing condition, K = Q, is one of the most widely used gauge conditions to fix the lapse 
function, a, during numerical evolution, since it has a feature of singularity avoidance. The condition will turn 
to an elliptic equation for a. However, in the case of Schwarzschild spacetime, we can express a maximally- sliced 
hypersurface at an arbitrary time without full numerical time integration. The recipe is given by Estabrook et al [30] , 
and we introduce the procedure briefly in Appendix B. 

By specifying a particular time t, where t is supposed to express the time coordinate on a maximally sliced 
hypersurface (see detail in Appendix B), we get metric components {a, (3 and 7ij). The procedure requires us to 
solve ODE, but it is not a result of time integration. We then calculate AFs as the previous ones. 

We show several snapshots in Fig. 7. We picked up 3-shces of f = 0, 1, 2, • • • , 5 for (a) the standard ADM and (b) 
Detweiler-type adjustment for the Schwarzschild coordinate. These initial data {i = 0) match those in Fig. 1(a) and 
Fig. 2(b), respectively. It is well known that the maximally-sliced hypersurfaces approach (or stop approaching to the 
singularity) rmin (3/2)M as f ^ oo. The snapshots here correspond to Vmin = 2.00, 1.90, 1.76, 1.66, 1.60 and 1.56 
in unit M. 

The figures show that AFs are changing along their evolution although not drastically, and also that their evolution 
behavior is dependent upon the choice of adjustments. The evolution makes the AFs' configuration converge, which 
is expected from the nature of maximal slicing. 

V. CONCLUDING REMARKS 

Motivated by performing a long-term stable and accurate numerical simulation of the Einstein equation, we pro- 
posed to adjust evolution equations by adding constraint terms by analyzing the constraint propagation equations in 
advance. The idea is to construct an asymptotically constrained evolution system, which is robust against violation 
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of the constraint. This method works even for the ADM formulation (which is not a hyperbohc system) against flat 
background spacetimc [25], and here we apphed the analyses to a curved spacetime, a spherically symmetric black 
hole spacetime. 

Recently, several numerical relativity groups report the effects of adjustments. They are mostly searching for 
a suitable combination of multipliers through trial and error. We hope our discussion here helps to understand 
the background mathematics systematically, though it may not be the perfect explanation. The main difference 
between our analysis and actual numerical studies is that this is a local analysis only on the evolution equations. The 
actual numerical results are obtained under a certain treatment of the boundary conditions, the choice of numerical 
integration schemes and grid structures, and also depend on the accuracy of the initial data. We think, however, 
that our proposal is an alternative to the hyperbolicity classification in that it includes the non-principal part of 
the evolution equations, and we expect that the discussion here will provide fundamental information on the stable 
formulation of the Einstein equations to the community. Although we have not shown any numerical confirmations 
of our conjecture in this article, we remark that the amplification factors explain the constraint violation behaviors 
quite well in the Maxwell equations and in the Ashtekar version of the Einstein equation on the flat background [22], 
and also that one of our examples explains a successful case of Kelly et al [11]. 

We presented a useful expression for analyzing ADM constraint propagation in general in Appendix A, and several 
analytic predictions of the adjustments for the Schwarzschild spacetime in Table. III. We searched when and where 
the negative real or non-zero imaginary eigenvalues of the constraint propagation matrix appear, and how they depend 
on the choice of coordinate system and adjustments. Our analysis includes the proposal of Detweiler (1987), which is 
still the best one though it has a growing mode of constraint violation near the horizon. 

We observed that the effects of adjustments depend on the choice of coordinate, gauge conditions, and also on its 
time evolution. Therefore our basic assumption of the constancy of the multipliers may be better to be replaced with 
more general procedures in our future treatment. We have already started to study this issue by applying the recent 
development of the computational techniques in classical multi-body dynamical systems [32], and hope that we can 
present some results soon elsewhere. 
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APPENDIX A: GENERAL EXPRESSIONS OF ADM CONSTRAINT PROPAGATION EQUATIONS 

For the reader's convenience, we express here the constraint propagation equations generally, considering the ad- 
justments to the evolution equations. We repeat the necessary equations again in order for this appendix to be read 
independently. 

1. The standard ADM equations and constraint propagations 

We start by analyzing the standard ADM system, that is, with evolution equations 

daij = -2aKij+Vi(3j+Vj(3i, (Al) 
dtKij = aRf) +aKKij-2aKikK''j-ViVja 

+{Vil3'')Kkj + {Vil3^)Ku + l3''VkKij, (A2) 

and constraint equations 

H := R^^^i + K'^ - KijK'\ (A3) 

Mi := VjKU-ViK, (A4) 

where {'jijjKij) are the induced three-metric and the extrinsic curvature, {a, Pi) are the lapse function and the shift 
covector, Vj is the covariant derivative adapted to ^ij, and r\^^ is the three- Ricci tensor. 

The constraint propagation equations, which are the time evolution equations of the Hamiltonian constraint (A3) 
and the momentum constraints (A4). 
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a. Expression using H and A4i The constraint propagation equations can be written as 

dtH = l3^{djH) + 2aKH-2a^'\diMj) 

+«(9,7„fc)(27'"S'^' - 7"S'^)>(j - ^l'^{dja)Mu (A5) 
dtMi = -{l/2)aid,n)-{d,a)n + f3^{djM,) 

+aKMi - P''-f'\daik)Mj + {diPkh^'Mj. (A6) 

This is a suitable form to discuss hypcrboUcity of the system. The simplest derivation of (A5) and (A6) is by using 
the Bianchi identity, which can be seen in Frittelli [26]. 
A shorter expression is available, e.g. 

= (3^Vin + 2aKn-2a{ViM^)-i{Via)M^, (A7) 
dtMi = -{l/2)a{d,n)-{d,a)n + (3^ViMi + aKMi + {V,dj)M^ 

= -{l/2)a{WiH) - {Wia)n + p^ViMi + aKMi + (ViA)A^', (A8) 

or by using Lie derivatives along an^, 

£an^n = 2aKH-2a^-^''^di{,/^M'-)-A{dia)M\ (A9) 
£ccn'-Mi = -{l/2)a{di'H)-{dia)n + aKMi. (AlO) 

h. Expression using ^ij and Kij In order to check the effects of the adjustments in (Al) and (A2) to constraint 
propagation, it is useful to re-express (A5) and (A6) using 7ij and Kij. By a straightforward calculation, we obtain 

an expression as 

dtH = H^{^an^n) + Hr"d,{dtJ,nn) + Hl''^''^,^J{^tJ„^n) + HridtK^n), (AH) 
dtMi = Mu""\dt^mn) + M2i'"''^dj{dt-/mn) + M^i"'^ {dtK^n) + Mu'^'^d.idtK^n), (A12) 



where 



-"1 • kj pi^ ^ T J- 1 

+7'^7"^(5^7'"')(a,7fcp) - 7'"V(5^7'^)(a,7fep) - 2i^if'"" + 2K^jK^\ (A13) 

-27™T" - (3/2)7*^ (5^7™") + 7'"^'(5j7™) + 7^^*, (A14) 

+ (Al5) 

2(iC7™" - i^""), (A16) 

^^^diK^j) - ^^^{djK^i) + {l/2){dj^^^)KU + T^K^i, (A17) 

-^rnjj^n. ^ {I I2)^rnn j^j . ^ (l/2)ii'™"5|, (A18) 

_^npm _ (1/2) (9^7™"), (A19) 

7'"-'(Jf -7'""(5^', (A20) 



Tj'iriin 

Tj-ijmn 
■"3 



where we expressed F™ = r-^7*-' . 



2. Adjustments 

Generally, we here write the adjustment terms to (Al) and (A2) using (A3) and (A4) by the following combinations, 

adjustment term of dtjij : +PijH + Q^ijMk +P^ij (Vfe7^)+g'='^,(Vfc^^^), (A21) 
adjustment term of dtKij : ^-RijU + S^ijMk + r^ij (VfeW) + s^\j {VkMi), (A22) 

where P,Q,R,S and p,q,r,s are multipliers (please do not confuse Rij with three Ricci curvature that we write as 
Rij^)- We adjust them only using up to the first derivatives in order to make the discussion simple. 
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By substituting the above adjustments into (All) and (A12), we can write the adjusted constraint propagation 

equations as 

dtH = (original terms) 

+Hr[Pmn'H + Q^muMk+P^mniVkn) + q"' mni^ kMl)] 
+Hr''di[Pmnn + QtnMk+P^mni^kU) + Q^^ mn{^ kMl)] 
+Hi^"'^didj[Pmnn + QL^fc +/m-(VfcW) + g'='™„(VfeA^O] 

+Hr[Rmnn + S'^^Mk + r'^mni^kH) + s''^ mn{^ kMl)], (A23) 

dt-M-i = (original terms) 

+ M2/""'d,[Prann + Oj;„Alfe + (VfeW) + g'^' ( VfcTW, )] 

+M3i'""[i^„^„W + Si^Mk + r'^mniykH) + s'='„^„(Vfc7W0] 

+Mj^"dj[Rmnn + Si^Mk + r'^mniVkH) + s'='„^„(VfeA4^)]. (A24) 

Here the "original terms" can be understood either as (A5) and (A6), or as (All) and (A12). Therefore, for example, 
wc can sec that adjustments to dt^ij do not always keep the constraint propagation equations in the first order form, 
due to their contribution in the third adjusted term in (A23). 

We note that these expressions of constraint propagation equations are equivalent when we include the cosmological 
constant and/or matter terms. 



APPENDIX B: MAXIMALLY SLICING A SCHWARZSCHILD BLACK HOLE 

The maximal slicing condition, ii' = 0, is one of the most widely used gauge conditions to fix the lapse function, 
a, during numerical evolution, since it has a feature of singularity avoidance. The condition will turn into an elliptic 
equation for a, such as 

V'Via = ( + K'^)a = K,jK'' a, (Bl) 

where we used the Hamiltonian constraint in the second equality. However, for the case of Schwarzschild spacetime, 
we can express a maximally-sliced hypersurface at an arbitrary time without full numerical integration. The recipe is 

given by Estabrook et al [30], and we introduce the procedure here briefly. 
The goal is to express Schwarzschild geometry with a metric 

ds^ = {-a^ + A-^p'^)dP + 2pdidr + Adr^ + r^{de^ + sin^ 9 d<fi^), (B2) 

where a(f, r), /3(t, r) and A{t,r), and we suppose i=const. on the maximally-sliced hypersurface. 
The original Schwarzschild time coordinate, t, is now the function of t{t,r), with 



dt 
dt 

dt 
dr 



= a\/A (B3) 

r=const. 

y/AT 1 



t=const. 1-2M/T 
where T is an arbitrary function of t, while A is given 



(B4) 



2M , 

using both Hamiltonian and momentum constraints. (B4) can be integrated as 

T f^^'^^ 1 1 

*=-tt/ ^ dx, (B6) 

M JM/r 1 - 2a; VI - 2a; + T'^M-^x^ ^ ' 

where X{T) is required to be a smaller real root of the quartic equation, 

1 - 2a; + T^M-V = 0, (B7) 
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in order to be consistent with (B3) under the boundary condition, 



the smoothness across the Einstein- Rosen bridge, r = r„ 

that turns dt/dr\r^r^^^ oo and T = T(r™„) = -\/''mj„(2M - rmin)- 
By identifying t t at spatial infinity, we get 

X{T) 



at i = 0, 



(B8) 



2a; v/Fr2^irY2M^^ 



dx, 



(B9) 



where the integration across the pole at x = 1/2 is taken in the sense of the principal value. 

In summary, if we specified a parameter T (0 < T < 3-\/3/4 corresponds to < f < oo), then we obtain the 
coordinate time t and f from (B6) and (B9) via (B7). Then we obtain the metric components by 



2M r2 
1 + -r 



1 + 



a{r,T) = 

l3{r,T) = aATr-'^, 
and (B5), where dT/dt in (BIO) can be calculated using 



Jo 



dx 



(1 - 2a; + T2M-4a;4)2/3 



(BIO) 

(Bll) 



dT 



dt 

— = lim 

dT Y^x 



1 



(1 - 2T2M-4X3)(2X - l)Vl-2y + T2M-4y4 
dx 



M Jo (1 - 2a; + T'^M-^x'^)^/^ 



(B12) 



They arc functions of r and the minimum value of r (that is, at a throat), rmin, is given by a larger real root of 



2Mr-^ + r2 = 0. Note that the fact r„ 



(3/2)Af when t ^ oo indicates the singularity avoidance property of 



maximal slicing condition. The comparison of this feature with the harmonic slicing condition is seen in [31]. 
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No. name 


adjustments (non-zero part) 


TRS 


1st? 


motivations 


standard ADM 


no adjustments 


- 


yes 




1 original ADM 


Rij = Kpajij 


yes 


yes 


kf = —1/4 makes original ADM 


2-a Detweiler 


(2.18)-(2.21), KL 


no 


no 


proposed by Detweiler [27] 


2-P Detweiler P-part 


Pij = —KLa^jij 


no 


no 


only use Pij term of Detweiler-type 


2-b Uetweiler h-part 




no 


yes 


only use S'" ij term of Detweiler-type 


2-s Detweiler s-part 




no 


no 


only use J^^j term of Detweiler-type 


3 simplified Detweiler 


Pij = -KLa-fij 


no 


no 


similar to above No. 2-P, but a monotonic 


4 constant ii', , 






yes 


us(h1 by I'ciiii Statt' gump ^ii^ 



TABLE I: List of adjustments we tested and plotted. (See more cases in Table. III). The column of adjustments are nonzero 
multipliers in terms of (2.13) and (2.14). The column 'TRS' indicates whether each adjusting term satisfies the time reversal 
symmetry or not based on the standard Schwarzschild coordinate. ('No' is the candidate that makes asymmetric amplification 
factors (AFs).) The column '1st?' indicates whether each adjusting term breaks the first-order feature of the standard constraint 
propagation equation, (2.11) and (2.12). ('Yes' keeps the system first-order, 'No' is the candidate of breaking hyperbolicity of 
constraint propagation.) 



adjustment \ coordinate 


standard (4.1) 


isotropic (4.2) 


iEF (4.3) 


PG (4.4) 




exact 


maximal 


exact 


tiEF = 


tpG=Q 


standard ADM 


Fig. 1(a) 


Fig.7(a) 


Fig.4(a) 


Fig.4(b) 


Fig.4(c) 


1 original ADM 


Fig. 1(b) 










2-a Detweiler 


Fig.2 


Fig.7(b) 


Fig.5(a) 


Fig.5(b) 


Fig.5(c) 


2-P Detweiler P-part 


Fig.3(a) 










2-S Detweiler S-part 


Fig.3(b) 










2-s Detweiler s-part 


Fig.3(c) 










3 simplified Detweiler 


Fig.3(d) 










4 constant Rrr 








Fig.6(a) 


Fig.6(b) 



TABLE II: List of figures presented in this article. The adjustments are explained in Table. I. 



16 



No. 


No. in 

Table. I 


adjustment 


1st? 


TRS 


Sch/iso coorc 
real. 


Is. 

imag. 


iEF/PG CO 
real. 


ords. 

imag. 










yes 












P-1 


2-P 


p. . — K r rv^^- ■ 


no 








Tnakps* 2 IMpff 


Tint annarpTif 


P-2 


3 


P- ■ — K- r ■ 


no 








Tnakpc; ^ ^pcr 


Tint annarpnt 


P-3 




P- - P — ic riv P — (c/v 






cl 1 o'n'I'l^T" on 1 l\l r>o' 
Oil tiiiiji y (^lii. 1 > '-'ti ■ 


IILJL CLUlJCll (.^lllj 




Tint SI T^nSI T"PTTt 


P-4 




P — K-V ■ 

-Tij ft, j'tj 


no 


no 






tnflkpQ 2 IMpff 


Tint annflrpnt 


P-5 




P — 

J ft, frr 






voH Pot; /pTil l\pcr 




t'pH PnQ /pTil IMpo" 


Tint PT^T^PTpnt 








no 




N/A 


N/A 


K ~ 1.35 min. vals. 


not apparent 


0-2 




n'' - ■ — K 

^ I J ^ rr — ri. 




yes 


red. abs vals. 


not apparent 


red. abs vals. 


not apparent 






(Qfc , . . — . or O"^ ■ ■ — AsYv^'' ■ 

'^t '* '^t 




vos 


red. abs vals. 


not apparent 


enl. Neg. 


enl. vals. 


O-i 




O'" , , 0' r, = K-, 


IICJ 


A'OS 


red. abs vals. 


not apparent 


red. alls vals. 


not apparent 


R-1 


1 


''J " f'J 


yes 


yes 


kf = —1/4 min. abs vals. 


kf = —1/4 min. vals. 


R-2 


4 




vcs 


no 


not apparent 


not apparent 


red. Pos. /enl. Neg. 


cnl. vals. 


R-3 




i-rr ' / rr 


vcs 


no 


cnl. vals. 


not apparent 


red. Pos. /enl. Neg. 


cnl. vals. 


S-i 


2-S 






ii( ) 


lint apparent 


not apparent 


not apparent 


iu)t apparent 


S-2 










makes 2 Neg. 


not apparent 


makes 2 Neg. 


not apparent 


n-1 




ir\^ ■ ■ Ti^ ■ ■ — icrv^- ■ 






red. Pos. 


red. vals. 


red. Pos. 


enl. vals. 


r)-2 




Ti^ ■ ■ Ti^ — /crv 
P zj p rr — rvLX 


no 




red. Pos. 


red. vals. 


red.Pos/enl.Neg. 


enl. vals. 


n-.3 

P o 




T>^ - - nrJ' *C/W 

P xj p rr — rvLX jrr 


no 




mates 2 Neg. 


enl. vals. 


red. Pos. vals. 


red. vals. 






kl rr , — K-rv-v- ■ 


no 


no 


K = 1/2 min. vals. 


red. vals. 


not apparent 


cnl. vals. 


a-2 




kl rr 

y y rr r\'(-i jrr 


no 


ves 


rod. abs vals. 


not apparent 


not apparent 


not apparent 


r-1 




k r 


no 


ves 


not apparent 


not apparent 


not apparent 


enl. vals. 


r-2 




^r flCX 


no 


yes 


red. abs vals. 


enl. vals. 


red. abs vals. 


enl. vals. 


r-3 




ij T" rr — txiOf^rr 


no 


yes 


red. abs vals. 


enl. vals. 


red. abs vals. 


enl. vals. 


s-1 


2-s 


8"^, - (l/3)7ij7''] 


no 


no 


makes 4 Neg. 


not apparent 


makes 4 Neg. 


not apparent 


s-2 




S ij S ij — l'\,Oi'^ij 


no 


no 


makes 2 Neg. 


red. vals. 


makes 2 Neg. 


red. vals. 


s-3 




S %j S ff — KOi'^ff 


no 


no 


makes 2 Neg. 


red. vals. 


makes 2 Neg. 


red. vals. 



TABLE III: List of adjustments we tested in the Schwarzscliild spacetime. The column of adjustments arc nonzero multipliers 
in terms of (2.13) and (2.14). The column '1st?' and 'TRS' are the same as Table. I. The effects to amplification factors (when 
K > 0) are commented for each coordinate system and for real/imaginary parts of AFs, respectively. The 'N/A' means that there 
is no effect due to the coordinate properties; 'not apparent' means the adjustment does not change the AFs effectively according 
to our conjecture; 'enl./red./min.' means enlarge/reduce/minimize, and 'Pos. /Neg.' means positive/negative, respectively. 
These judgements are made at the r ~ O(IOM) region on their t = slice. 
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(a) 

no adjustments (standard ADM) 




5 10 15 20 5 10 15 20 



FIG. 1: Amplification factors (AFs, eigenvalues of homogenized constraint propagation equations) are shown for the standard 
Schwarzschild coordinate, with (a) no adjustments, i.e., standard ADM, (b) original ADM {kf = —1/4) and (c) an adjusted 
version with different signature, kf = +1/4 [see eq. (2.17)]. The solid lines and the dotted lines with circles are real parts and 
imaginary parts, respectively. They are four lines each, but actually the two eigenvalues are zero for all cases. Plotting range 
is 2 < r < 20 using Schwarzschild radial coordinate. We set fc = 1,/ = 2, and m = 2 throughout the article. 




FIG. 2: Amplification factors of the standard Schwarzschild coordinate, with Detweiler type adjustments, (2.18)-(2.21). 
Multipliers used in the plot are (a) kl ~ +1/4, (b) kl ~ +1/2, and (c) kl = —1/2. Plotting details are the same as Fig.l. 
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FIG. 3: Amplification factors of the standard Schwarzschild coordinate, with Detweiler type (partial contributions) adjust- 
ments. Fig. (a) is the case No. 2-P in Table. I, adjusting with Pij = —KLa'^'fij , i.e. (2.18), and else zero. Similarly, (b) is No. 
2-S, S^ij-part (2.20) and else zero, and (c) is No. 2-s, s'°'ij-part (2.21) and else zero. Fig. (d) is No. 3, Pij — —KLdJij and 
else zero, which is a minor modification from (a). We used kl = +1/2 for all plots. Plotting details are the same as Fig.l. 
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(a) 

isotropic coordinate, no adjustments (standard ADM) 




FIG. 4: Comparison of amplification factors between different coordinate expressions for the standard ADM formulation (i.e. 
no adjustments). Fig. (a) is for the isotropic coordinate (4.2), and the plotting range is 1/2 < Tiso- Fig. (b) and (c) are for the 
iEF coordinate (4.3) and the PG coordinate (4.4), respectively, and we plot lines on the t = Q slice for each expression. [See 
Fig. 1(a) for the standard Schwarzschild coordinate.] The solid four lines and the dotted four lines with circles are real parts 
and imaginary parts, respectively. 
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FIG. 5: Similar comparison with Fig. 4, but for Detweiler adjustments, kl = +1/2 for all plots. [See Fig. 2(b) for the standard 
Schwarzschild coordinate.] 
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FIG. 6: Amplification factors of the adjustment No. 4 in Table. I. = +1/2, -1/2 for iEF/PG coordinates. [See Fig.4(b)(c) 
for the standard ADM system.] 
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(a1) (a2) 

maximally sliced evolution, standard ADM maximally sliced evolution, standard ADM 




FIG. 7: Amplification factors of snapshots of maximally-sliced evolving Schwarzschild spacetime. Fig (al) and (a2) are of the 
standard ADM formulation (real and imaginary parts, respectively), (bl) and (b2) are Detweiler's adjustment {kl = +1/2). 
Lines in (al) and (bl) are the largest (positive) AF on each time slice, while lines in (a2) and (b2) are the maximum imaginary 
part of AF on each time slice. The time label t in plots is t in Appendix B. The lines start at Vmin — 2 {t = 0) and rmin ~ 1.55 
(t = 5). 



